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Abstract 

Recent developments in quantum technology mean that is it now possible to manipulate systems and 
measure fermion fields (e.g. reservoirs of electrons) at the quantum level. This progress has motivated 
some recent work on filtering theory for quantum systems driven by fermion fields by Korotkov, Milburn 
and others. The purpose of this paper is to develop fermion filtering theory using the fermion quantum 
stochastic calculus. We explain that this approach has close connections to the classical filtering theory 
that is a fundamental part of the systems and control theory that has developed over the past 50 years. 

1 Introduction 

A basic problem in control and communication systems is that of extracting information from a signal that 
may contain noise. Problems of this kind are known as filtering problems. One common scenario concerns 
the problem of estimating variables of a system from partial, and typically noisy, information. Here, the 
word 'system' refers to the entity of interest, which may be a machine being controlled, or it may be a signal 
model. As remarked in [1] sec. 1.2], filtering theory developed in response to demands from applications. 
For example, the Kalman filter [22 was developed at a time when significant efforts were underway in 
aerospace engineering in the early 1960 's. The Kalman filter is the solution of a filtering problem based on 
a statistical model involving Gaussian stochastic processes and linear dynamics (statistical filtering dates 
back to Kolmogorov [33J and Wiener [35]). A more general theory of nonlinear filtering was developed 
during the 1960's: Kushner [27], Stratonovich [3S], Duncan [ID], Mortensen [32"] . and Zakai [33J. These 
filtering results boil down to determination of conditional expectations in dynamical contexts, and this 
may be achieved using powerful tools from stochastic calculus, including idealized Wiener process models 
for noise and Ito stochastic differential equations, [HI Chapter 18]. 

At the present time, developments in quantum technology are demanding methods for statistical esti- 
mation (among other things). Quantum technologies are those technologies that depend on the laws of 
quantum mechanics for their operation. Examples of quantum technologies currently under development 
include quantum computers and atom lasers. Progress in quantum filtering has been slow, due to both 
conceptual and practical experimental issues concerning the measurement of quantum systems. However, 
significant advances were made in the field of quantum optics, a field of study concerned with quantum 
properties of light and the interaction of light and matter. In quantum optics, thanks to the invention 
of the laser and other experimental developments, quantum effects became accessible and a sophisticated 
theory of open systems, measurement theory, and quantum noise emerged. In particular, we mention the 

'Institute for Mathematics and Physics, Aberystwyth University, SY23 3BZ, Wales, United Kingdom. Email: 
jug@aber.ac.uk 

^School of Mathematical Sciences, University of Nottingham, Nottingham, NG7 2RD, United Kingdom. Email: 
madalin . gut a@ not t ingham .ac.uk 

* School of Engineering, Australian National University, Canberra, ACT 0200, Australia. Email: 
Matthew.James@anu.edu.au 

^School of Engineering, Australian National University, Canberra, ACT 0200, Australia. Email: Hen- 
dra. Nurdin@anu.edu . au 



1 



quantum filtering results due to Belavkin [4] , [5] that extend the thinking behind the above-mentioned clas- 
sical filtering theory (Belavkin was a student of Stratonovich in the 1960's). Belavkin employed quantum 
stochastic methods that were developed in the 1980's to describe quantum noise and a quantum generaliza- 
tion of the Ito calculus, (20] , [15] . Belavkin's results involved an implementation of conditional expectation 
in a dynamical quantum mechanical context, [7J. We also mention related work in quantum optics by 
Carmichael [8] and Wiseman and Milburn [41], [40] which used more direct methods, and employed the 
terms 'stochastic master equation' and 'quantum trajectories' in connection with quantum filtering. 



atom detector 




Figure 1: Schematic representation of the detection of light emitted by an atom. The measurement signal 
Y(t) (integrated photocurrent) is proportional to the number of photons in the light incident on the 
detector up to time t. 

The absorption and emission of light by atoms was one of the earliest problems studied in quantum 
optics, [TT], J5S] Chapter 1], Chapter 1]. Light is a type of electromagnetic field, which in quantum 
mechanics is described as a boson field, with characteristic Bose-Einstein statistics. Photons are the well- 
known 'particles' or 'quanta' of light - an example of a boson. Boson fields may be used in other contexts, 
such as in the description of vibrations in solid materials, or in the modeling of dissipation to a heat bath. 
To give some indication of how quantum filtering may be used in quantum optics, consider the setup of 
Figure [l] which illustrates the detection of light emitted by an atom. If we model the atom as a two- level 
system with ground and excited states, then the occupation number operator n for the excited state is a 
physical observable of interest (a self-adjoint operator with eigenvalues and 1). The event corresponding 
to the eigenvalue 1 means that the atom is in the excited state and so contains one 'quanta' of energy. 
Any energy in the atom may be emitted to the ambient field, a dissipative process. Indeed, the mean 
value n(t) = e _7 *n(0) decays exponentially (7 is a parameter describing the strength of the coupling of 
the atom to the light field). The photodetector (PD) shown in Figure [l] is taken to be an idealized device 
that produces a classical (i.e. non-quantum) photocurrent that is proportional to the number of photons 
in the field and may be processed using conventional analog or digital electronics. A quantum filter is such 
a processing system that is designed to provide estimates of, in this case, atomic observables (which are 
not directly accessible). The quantum filter for the conditional expectation h(t) of the occupation number 
given the photocurrent is 

dn{t) = -~fh(t)dt - h(t)(dY(t) - jri(t)dt), (1) 

where Y(t) is the integrated photocurrent signal. This filter has a form that is familiar to control engineers, 
which combines a prediction term with an update term involving the innovations process W(t) defined 
by dW(t) = dY(t) — ^h(t)dt. We refer the reader to the book [35] and the tutorial paper [7] for further 
information on quantum filtering involving boson fields. 

In quantum field theory, there is another type of field that is distinct from boson fields. These are called 
fermion fields, and the quanta of these fields are called fermions, the electron being an important example. 
Fermion fields have Fermi-Dirac statistics. Fermion fields arise in mesoscopic systems, such as quantum 
dots, which are of considerable technological importance (e.g. for use in quantum computers). Quantum 
dots may be fabricated in semiconductor materials to confine one or a few electrons to a region the size 
of a few nanometers. Fermion fields may be used to describe the flow of electrons at the quantum level. 
Figure [2] provides a schematic representation of a quantum dot connected to two fermion field channels 
representing a source and a sink. Experimentally it has been much harder to extract quantum properties 
of mesoscopic systems compared to quantum optical systems, however, it is now experimentally feasible 
to do so. 
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Figure 2: Schematic representation of quantum dot, through which current tunnels from the source (L) 
ohmic contact to the drain (R), [31]. 



Some important results have been obtained concerning quantum filtering theory for the case of fermion 
fields. Korotokov [23], [25] developed a phenomenological approach using Bayesian methods, and Goan 
et al [T7], [TB] adapted methods from quantum optics. The purpose of the present paper is to develop 
quantum filtering theory for systems driven by fermion fields using the fermion quantum stochastic calculus, 
Applebaum and Hudson, [2], Hudson and Parthasarathy [2T], Milburn [3T], Gardiner [T3]. We explain 
that quantum filtering, involving both boson and fermion fields, and classical filtering share a common 
mathematical foundation in conditional expectations, and that stochastic calculus provides powerful tools 
with which quantum and classical filtering problems can be solved. The key difference between the classical 
and quantum problems is that in the quantum cases, the random variables and stochastic processes (such as 
quantum observables and field operators) have non-commutative algebraic structures that are fundamental 
to quantum mechanics (classical variables are represented as scalar valued functions and so commute under 
pointwise multiplication). The principle algebraic distinctions between the boson and fermion cases are 
the commutation and anticommutation relations satisfied by these fields (which underly the Bose-Einstein 
and Fcrmi-Dirac statistics), and the need to take parity into account in the fermion case. 

In an effort to be as concrete as possible, in this paper we develop the filtering theory using a model 
with just enough generality to solve filtering problems for two examples. The first example we consider 
(see Section 5.1 ) concerns the quantum dot shown in Figure [2j where the electron flow in the right channel 
is monitored, and this information is used to estimate the number of electrons (0 or 1) in the dot. Our 
second example is based on a more detailed model for the process of photodetection shown in Figure [T] 
[161 sec. 8.5]. This more detailed model is shown in Figure |3j and involves one boson and two fermion 
field channels. 
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Figure 3: Schematic representation of a model for the detection of the photons emitted by an atom. 
The model includes one boson field channel B(t), two fermion field channels Ao(t) and Ai(t), a two- level 
system to describe the atom, and a three-level system capturing the essential behavior of the detector. 
The notation is defined in Section 15.21 

The paper is organized as follows. In Section [2] we review some basic ideas about classical filtering 
theory, and in particular, we summarize how the fundamental nonlinear filtering results may be obtained 
using classical stochastic calculus. Section [3] provides a description of the quantum stochastic model 
that is needed to formulate the filtering problem. This section include a brief review of some aspects of 
quantum mechanics, with an emphasis on describing the boson and fermion commutation relations. The 
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main fermion filtering results are presented in Section |4j and the examples are given in Section [5j Two 
appendices briefly summarize some basic aspects of classical and quantum stochastic calculus, and show 
how parity arises in systems with fermionic degrees of freedom. 

Notation: In what follows the symbols E and E represent classical and quantum expectations, respec- 
tively. The commutator of two operators A and B is denoted [A, B] = AB — BA, while the anticommutator 
is written as {A,B} = AB + BA. The Dirac notation for a vector ip in a Hilbert space Sj is and 
the inner product is written as ((j),tp} or (<j>\ip). The quantum expectation for an observable X when the 
system is in a state described by the vector ip is E</,[X] = {ip, Xip) or (tp\X\ip). The adjoint of an operator 
X is denoted by X* . 



2 Classical Filtering Theory 

In this section we review some of the fundamental concepts and results concerning classical (non-quantum) 
filtering theory that will assist with understanding the quantum filtering results to be presented below. Let 
£ and Y be classical random variables, with joint density p^ t y(x,y)- In the absence of any measurement 
data, one may simply use the marginal density p^(x) = J p^ : Y(x,y)dy to make inferences about £; for 
instance, the mean £ = J p^(x)dx gives us an indication of the value an observation of £ may yield. If a 
value y of Y is observed, then the density for £ is revised to the conditional density 

reflecting an increase in knowledge. The conditional mean is defined to be 

i= / xp ilY (x\y)dx, (3) 



which we note is a function of the observed data y. More generally, if <fi is an arbitrary function, we may 
compute the conditional expected value 

$ = J <t>{x)px\y{x\y)dx (4) 

of the random variable The RHS of Q is an explicit expression for the conditional expectation, 

which is denoted more generally as 

0=7r(0)=E[0(oin (5) 

As a simple example, suppose that 

Y = £ + V, (6) 

where £ and V are independent Gaussian random variables with means £ and 0, and variances and 
1. The expression Q for the observations Y is an instance of the fundamentally important "signal plus 
noise" models widely employed in control and communications systems. The conditional mean £ is also 
Gaussian p] and is given by 

£ = t + Xt(l + Zt)- 1 (Y-£). (7) 

This expression shows how the conditional mean updates the prior mean £ by the addition of a term 
W = Y — £, called the innovation. The innovation represents the new information about £ that is gained 
from an observation of Y. 

The innovation is related to the minimum variance property of the conditional mean, which means 
that £ minimizes the variance of the error "error" £' = £ — £ over all estimators £. Here, estimators £ are 
random variables that arc functions of the observation Y, that is, random variables that belong to the 
subspace generated by Y^ The least squares property has a nice geometrical interpretation, where £ is 
the orthogonal projection of £ onto W, Figure |4j 

lr rhat is, 'Sf is the subspace of square integrable random variables that are measurable with respect to the cr-algebra o{Y) 
generated by Y . 
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Figure 4: The conditional expectation £ = E[£|Y] is the orthogonal projection of X onto the subspace 
[D Fig. 5.2-1]. 



In general, if <f> and Y have well defined expectance then the conditional expectation E[^|Y] is defined 
to be the unique (up to a set of measure zero with respect to the underlying probability measure) random 
variable 4> = E[<t>\Y] € such that [44, Chapter 1] 

E[#] = E[#] (8) 

for all bounded random variables 4> G *3f . 



2.1 Nonlinear Filtering of Classical Systems in Continuous Time 

In control and communications systems, signals are often modeled as stochastic processes, which are 
sequences of random variables. For linear Gaussian systems, the Kalman filter computes the conditional 
mean and covariances in a causal manner as time progresses, as described in the book [1] for discrete time 
systems. Of interest to us here are general approaches to filtering in continuous time that exploit the 



power of the stochastic calculus, [T2J Chapter 18]. We refer the reader to Appendix A.l for some basic 
concepts concerning stochastic integrals and the Ito rule. 

Consider the stochastic system expressed as a system of Ito stochastic differential equations 

dm = giaWt + dV^t) (9) 
dY{t) = h(£(t))dt + dV 2 (t) (10) 

where Vi(t) and V<i(i) are independent Wiener processes. Here, £(t) represents a signal of interest that is 
not directly accessible to observation. Instead, a signal Y(t) is observed. 

Let us first consider the dynamics of expected values for this system. For any sufficiently regular 
function <fi write jt(4>) = <fi(£(t)) and define 

= E\j t (<t>)] (11) 
for the mean of the random variable j t (0) . Now by the Ito rule we have 

Mfi =^3)dt + j t (^)dV 1 (t))+j t (^)dt, (12) 



and so 



where 



|/*t(0 = (13) 

„ . , N 1 d 2 d> dd> ... 
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is the generator of the Markov process given by The Kolmo gorov equation for the density p(x, t) 
defined by ^t(4>) = J <fi(x)p(x,t)dx is 

J^ = £»> ( 15 ) 

where C*(p) = \£,p - ^(gp). 

Now suppose we wish to determine the differential equation for the conditional expectation 

n t (if>)=E[i t (4>)\SQ, (16) 

where Wt is generated by the observations Y(s), < s < t. For convenience, we write h t for h(£(t)). 



We might expect the filter equation to be a modification of the mean equation (13) along the lines of 
([7]). There are several commonly used methods for finding the filter equations, including the martingale 
approach (which lead to the Kushner-Stratonovich equation [27], |36j). the reference probability method 
(giving the Duncan-Mortensen-Zakai equation [ID] . [33] , [35]). and the characteristic function method, [35], 
[BJ, [3] which we will use in this paper. 

We suppose that the filter has the form 

d-K t {<t>) = F t (<t>)dt + G t (<j>)dY(t) (17) 

and consider, for any square- integrable function /, the stochastic process Ct(f) defined to be the solution 
of 

dC t (f) = f(t)C t (f)dY(t), C f (0)=I. (18) 
By the definition of conditional expectation (recall ([8])), we have 

m^)c t {f)] = n^M)CtU)] (19) 

for all /. We note that the Ito product rule implies I + II + III = where 

/ = E[{dK t y>)-4j t (4>))C t {f)], 

II = E [(7r t (4>) — j t (4>)) dC t (j)] , 
III = E [(dir t (<f) — dj t (<fi)) dCt(f)] ■ 

For the model above, we then have 
I = E [{F t (</>) + h t G t (0) - j t (£0)} C t (/)] dt = E [{F t (4>) + Tr t (h t )G t (0) - ir t (£0)} C t (/)] dt, 

II = E [{7T t (0) - j t (<f>)} / (t) C* (/) /l t ] tft = / (t) E [{7T t (<£) tt 4 (/i t ) - j t (# { )} C t (/)] tfc, 

JIJ = E[G t ((f>)f(t)Ct(f)}dt. 

As /(i) was arbitrary, we may separate the / independent and dependent terms to obtain 

F t {4>) = TTt (C<j>) - Tit (h t ) G t , G t ((/>) = TTt (<j>h t ) - TT t (<t>) TTt (ht) 

so that 

dTTt (0) = TTt (C(j)) dt + {TTt (#t) " 7T t (0) TTt {ht)} dW t 

where the innovations process is a ^ martingale (actually a standard Wiener process) given by 

dW(i) = dY(t) - Tr t (h t )dt. (20) 

The conditional density p(x, t) defined by / <fi(x)p(x,t)dx = E[j t (^)|^] satisfies the equation 

dp = C*(p)dt + (h-TT t (h))pdW(t). (21) 
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In the special case of linear systems, say <?(£) = a£ and h(£) = c£, with initial Gaussian states, the 
process £(t) is Gaussian, with mean = E[£(t)] satisfying the equation 

|(t) - af(t), (22) 

and variance r(t) = E[(£(t) — £(i)) 2 ] satisfying 

f (t) = 2aT(t) + 1. (23) 

The conditional mean = E[£(t)|$t] is Gaussian and is given by the Kalman filter equations 

d£(t) = a£(t)dt + cE(i)(dY(t)-c£(t)dt) (24) 
t(t) = 2a£(i) + 1 + c 2 £ 2 . (25) 

Here, S(t) = E[(f(t) - £(i)) 2 |^] is the conditional variance, a deterministic quantity (a special feature of 
the linear-Gaussian case). 



3 Quantum Stochastic Model 

In this paper we are interested in a quantum system interacting with quantum fields, for instance as 
sketched in Figure [5j Here the "box" represents a quantum system with finitely many degrees of freedom, 
such as an atom or a quantum dot or a photodetector. The input/output lines are quantum fields, 
representing reservoirs of electrons or photons coupled to the system. The model for this system has a 
natural input-output structure, with an input being the incident part of the field, while the output is the 
reflected part of the field which carries away information about the system. Our main goal in this paper 
concerns estimation of system variables given the results of monitoring the output of fermion channel 0. 
The purpose of this section is to describe a quantum mechanical model for this system. 
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Figure 5: Schematic representation of a system coupled to boson B and fermion A , A\ fields. 



3.1 Quantum Mechanics 

Quantum mechanics was developed in the 20th century in order to explain the behavior of light and 
matter on a small scale, |29j . [5]. Central to quantum mechanics are the notions of observables X, which 
are mathematical representations of physical quantities that can (in principle) be measured, and statej^] 
vectors ip, which summarize the status of physical systems and permit the calculation of expected values 
of observables. State vectors may be described mathematically as elements of a Hilbert space ft, while 
observables are self-adjoint operators on f). The expected value of an observable X when in state ip is 
given by the inner product (tp, Xip). 

2 The word 'state' is heavily overloaded in the physical and engineering sciences, though its core meaning as a way of 
minimally storing dynamical and statistical information is common to all interpretations. 
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The simplest non-trivial quantum system has two energy levels and is often used to model ground 
and excited states of atoms. Since the advent of quantum computing, this system is also known as the 
qubit, the unit of quantum information. The Hilbert space for this system is fj = C 2 , the two-dimensional 
complex vector space. The space of all operators is spanned by the Pauli matrices HH sec. 2.1.3], [THl sec. 
9.1.1]: 

, ( 1 \ / 1 \ / -i \ ( 1 



v o i J ' x y i o y ' u \ i o / ' z \ o -l 

Since in general operators do not commute, the commutator 

[A, B] = AB - BA, (26) 

which is a measure of the failure to commute, is frequently used. The basic commutation relations for the 
Pauli matrices are [o~ x ,a y ] = 2ia z , \a y ,a z \ = 2ia x , and [a z ,a x ] = 2ia y . 

The energy levels correspond to the eigenvalues ±1 of a z , and the corresponding eigenstates are referred 
to as the ground | — 1) and excited states |1). These eigenstates form a basis for F). In quantum mechanics, 
state vectors tp are normalized to one, so that we may write 

= «| - 1) +/3|1>, (27) 

where \a\ 2 + \(3\ 2 = 1. The operators a± = \{p x ± icr y ) are known as raising and lowering operators, 
and have actions er+| — 1) = |1) and <J-|1) = | — 1), as illustrated in Figure |6] (a). Thus <r + corresponds 
to the creation of a quanta of energy in the system, while <7_ corresponds to destruction. The number 
operator n = (7 + er_ counts the number of quanta in the system, which is in this case is or 1 (ground 
or excited). Note that a+ is the adjoint of er_: a+ = a*_. The raising and lowering operators satisfy the 
anticommutation relation 

{*-,<,*_} = I, (28) 

where the anticommutator is defined by 

{A,B}=AB + BA. (29) 

Note also that a 2 _ = 0. 

The postulates of quantum mechanics state that if an observable A is measured, the allowable mea- 
surement values are the eigenvalues of A. If the system is in state ip, the probability of observing the 
outcome Xj is 

Prob(A J ) = (il>,Pji>), (30) 

where Pj is the projection associated with the eigenvalue Xj. When an eigenvalue A-,- is recorded, the state 
"collapses" to 

^ ^/Prob(aj) ^ ^ 

This brief discussion of quantum measurement theory has focused on ideal measurements when the 
system is in a pure state (a state vector). Density operators p (self-adjoint non- negative operators of trace 
one) provide a more general notion of state, where the probability of outcome is given by Prob(Aj) = 
trfpPj], and the collapsed state is p' = p^^^ry- A theory of generalized measurements was developed 
largely during the 70's an 80's that allows for non-ideal circumstances [531 Chapters 2 and 8], and by 
the addition of ancilla systems all of quantum measurement theory can be seen to be a consequence of 
quantum conditional expectation, [7|. 

Another fundamental postulate of quantum mechanics states that the dynamics of states and observ- 
ables is unitary. That is, state vectors and observables evolve according to 

il>(t) = U(t)4>, X(t) = U*(t)XU(t), (32) 
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where the operator U(t) is unitary (U*(t)U(t) = U{t)U*(t) = I). The unitary U(t) is the solution of the 
Schrodinger equation 

jU(t) = -iHU(t), (33) 



where the observable H is called the Hamiltonian. The expressions in ( 32 1 provide two equivalent descrip- 
tions (dual), the former is referred to as the Schrodinger picture, while the latter is the Heisenberg picture. 
The differential equations in the Schrodinger and Heisenberg pictures are respectively 

j t m = -mm (34) 



and 



where H(t) = U*(t)HU(t). 



jX{t) = ~t[X(t),H{t)], (35) 



• \ t 

• I t 

- I t 

t - ■ I t 



(») (b) 

Figure 6: Energy level diagrams, (a) Two-level atom (qbit). (b) Harmonic oscillator. 

Another basic example is that of a particle moving in a potential well, [|29, Chapter 14] . The position 
and momentum of the particle are represented by observables Q and P, respectively, defined by 

(Qi>)(q) = #(<z), (PV0(<Z) = -*^(g) (36) 

for ip G — _L 2 (R). Here, q €E R represents position values. The position and momentum operators satisfy 
the commutation relation [Q, P] — i. The dynamics of the particle is determined by the Hamiltonian 
H = y — h ^muj 2 Q 2 (here, m is the mass of the particle, and ui is the frequency of oscillation). 

Energy eigenvectors ip n are defined by the equation Hijj n = E n ip n for real numbers E n . The system 

has a discrete energy spectrum E n = (n + n = 0, 1, 2, The state ipo corresponding to Eq is called 

the ground state. The annihilation operator 



wa>,_ P „ ,„„. 
2 2mw 

and the creation operator a* lower and raise energy levels, respectively: aip n = \fnif) n —\, and a*ip n = 
\/n + lip n +i, see Figure [6] (b) . They satisfy the canonical commutation relation 

[a,a*] = l. (38) 

In terms of these operators, the Hamiltonian can be expressed as H — uo(a*a + 4). Here, i 



a a 



is the number or counting operator, with eigenvalues 0,1,2,.... Using (35), the annihilation operator 



evolves according to 4;a(t) = — iwa(t) with solution a(t) = e~ lut a. Note that also a*(t) — e lut a*, and so 
commutation relations are preserved by the unitary dynamics: [a(t),a*(t)] — [a, a*] — 1. Because of the 
oscillatory nature of the dynamics, this system is often referred to as the quantum harmonic oscillator. 
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3.2 Boson and Fermion Fields 



We may think of observables as quantum random variables, and the key distinction with classical prob- 
ability is that quantum random variables do not in general commute. Indeed, if (fi, P) is a classical 
probability space then classical bounded real- valued random variables in L°°(Q, IF, P) have an interpreta- 
tion as multiplication operators that map the Hilbert space L 2 (Vl,F,P) to itself. Since all such operators 
commute with one another, bounded classical real-valued random variables are thus isomorphic to (and 
can be viewed as) commuting observables on L 2 (il, T , P); see [7j for further discussions, including the case 
of unbounded classical random variables. 

In this section we turn to the notion of quantum stochastic processes which are used to provide tractable 
models for how quantum systems evolve in dissipative environments, such as a quantum dot in an electron 
field, or an atom in an electromagnetic field. A detailed introduction to quantum fields is beyond the scope 
of the present paper, and we refer interested readers to the references for more information, particularly 
[29l Chapter 21], [35l Chapter II], Section 4]. Here we give some basic ideas needed for the filtering 
results to follow. 

In quantum field theory, a one dimensional quantum field (with parameter t) consists of a collection 
of systems each with annihilation a(t) and creation operators a*(t) used to describe the annihilation and 
creation of quanta or particles at index location or point t. a(t) and a*(t) are referred to as field operators, 
the annihilation and creation field operators, respectively. The index t may represent a range of variables, 
including position, frequency and time, and we assume here that t lies in a continuous interval T in R. 
Basic considerations lead to the postulate that the annihilation and creation operators must satisfy either 
the commutation relations 

[a(t),a*(t')] = 6(t-t'), (39) 

or the anticommutation relations 

{a(t),a*(t')} = S(t-t'), (40) 
for all t, t' € T, where S(t) denotes the Dirac delta distribution. 



Fields that satisfy the commutation relations (391 are called boson fields (e.g. photons), while fields 



that satisfy the anticommutation relations (40) are called fermion fields (e.g. electrons). In this paper we 
will take the parameter t to be time and T = [0, oo). In this case a(t) has the interpretation of annihilation 
of a photon (in the case of a bosonic field) or electron (in the case of fermionic field) at time t, whereas 
a* (t) has the interpretation of creation a photon (in the case of bosons) or electron (in the case of fermions) 
at time t. One can imagine these fields as a continuous collection or stream of distinct quantum systems 
(one quantum system for each t) hence, informally, quantum fields can be defined on some continuous 
tensor product Hilbert space % — ®te[Q,co)H-t, where T-L t is a Hilbert space for each t (of the quantum 
system arriving at time t). Although such an object can be rigorously defined and constructed, from 
a mathematical viewpoint it is much easier not to work directly with the field operators a(t) and a*{t) 
but with their integrated versions, the so-called smeared quantum field operators, as will be discussed 
below. Smeared quantum field operators can be constructed on Hilbert spaces known as Fock spaces 
(symmetric Fock space 'Ssym for bosons and antisymmetric Fock space $ a ntisym for Fermions) which have 
the character of a continuous tensor product Hilbert space. Modulo the specification of the statistics of the 
field, a quantum field has the character of a quantum version of white noise, while its integrated version 
can be viewed as a quantum independent increment process. Thus, exploiting the properties of smeared 
quantum fields, Hudson and Parthasarathy |20j were able to develop a quantum stochastic calculus which 
is essentially a quantum version of the Ito stochastic calculus. 

The model we use to describe the system shown in Figure [5] employs boson and fermion fields b(t) 
and a(t), respectively, parametrized by time t € [0, oo) which accounts for the time evolution of fields 
interacting with the system (e.g. an atom or quantum dot) at a fixed spatial location. In the remainder 
of this section we describe the quantum stochastic calculus that has been developed to facilitate modeling 
and calculations involving these fields, [20], [2], [15], [21], [35], [16], [7], [42]. Some basic aspects of quantum 
stochastic integrals and the quantum Ito rule are discussed in Appendix | A. 2 1 

The boson field channel B in Figure [5] is defined on a symmetric Fock space $ S ym- The commutation 
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relations for the boson field are [b(t), b*(t')} = 8(t — t'), from (39). For a boson channel in a Gaussian state, 
the following singular expectations may be assumed: 

{b*(t)b(t'))=N8(t-t'), {b{t)b*(t')) = (N + l)6(t-t'), (41) 
(b{t)b(t')) = M5(t - t'), (b*(t)b*(t')) = M*S(t - t'). (42) 

Here (X) is a standard notation used to denote the quantum expectation of a system operator X 29, 3 
(i.e., (X) = ELY]), TV > is the average number of bosons, while M describes the amount of squeezing in 
the field state. We have the identity |M| 2 < N(l + N). For a thermal state, M = and 

N = eP (E-\ ) _ V (43) 

where f3 = is the inverse temperature, E is the energy, and fi is the chemical potential. 

In this paper we will assume N = M = 0, which corresponds to the case of a boson field in the vacuum 
(ground) state. The vacuum boson field is a natural quantum extension of white noise, and may be 
described using the quantum Ito calculus. In this calculus, the integrated field processes B(t) — f Q b(s)ds 
(annihilation), B*(t) = f*b*(s)ds (creation) and A(t) — f*b*(s)b(s)ds (counting) are used. The non-zero 
Ito products for the vacuum boson field are 

dA{t)dA(t) = dA(t), dA{t)dB*{t) = dB*(t), dB(t)dA (t) = dB(t), dB(t)dB*(t) = dt. (44) 

We now specify the fermion channels Aq and A\ in Figure [5] We assume the followings singular 
expectations for a fermion field A, defined on an antisymmetric Fock space 'Santisym' 

(a*(t)a(t')) = NS(t-t'), (a(t)a*(t')) = (1 - N)S(t - t'), (45) 
(a(t)a(t')) = MS(t-t'), (a*(t)a*(t')) = M*5{t-t'). (46) 

In general we have < N < 1 along with the identity \M\ 2 < N(l - N). For a thermal state we have 
M = 0, and 

In what follows we take the zero temperature limit T — > 0. For fermion channel 1 we assume the energy 
is such that E < \l and so in the zero temperature limit this channel is fully occupied, N — 1, and the Ito 
rule 

dAl(t)dA!{t) = dt (48) 

applies for the corresponding integrated processes Ai(t) = f*ai(s)ds, A\(t) = J*al(s)ds. For fermion 
channel we fix E > /i, in which case N = 0, describing a reservoir which is unoccupied. The number 
process Ao(t) — J Q ao(s)ao(s)ds is well defined for fermion channel (but not for channel 1), and the Ito 
table is 

dA (t)dA (t) = dA (t), dA Q (t)dA* (t) = dA* Q {t), dA (t)dA {t) = dA {t), dA {t)dA* Q {i) = dt. (49) 
The fermion channels are defined on distinct antisymmetric Fock spaces Santisym' Santisym- 

3.3 System Coupled to Boson and Fermion Fields 

The system S illustrated in Figure [5] is defined on the Hilbert space Sjs, and so the complete system 
coupled to the boson and fermion fields is defined on the tensor product Hilbert space 

= $)S ® $sym ® fiintisym ® ^Lntisym- (50) 

Due to the presence of fermion field channels, it is necessary to introduce a parity structure on the 
collection of operators on this tensor product space, as explained in Appendix iBj We therefore have a 
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parity operator r on such that for all operators X and Y on we have t(XY) = t{X)t(Y) and 
t(X*) — t(X)* . Operators X such that t(X) = X are called even, while those for which t(X) — —X are 
called odd. Fermion annihilation and creation operators are odd, while the fermion number operator is 
even. All boson operators are even. A system operator, i.e. an operator X acting nontrivially on ios only, 
that is even will commute with all field operators, while an odd system operator will anticommute with 
odd fermion field operators. All boson field operators commute with all system operators and all fermion 
field operators. 

The Schrodinger equation for the complete system is 

dU(t) = {(S-I)dA(t) + dB*(t)L-L*SdB(t)--L*Ldt 

+dA* 1 (t)L 1 - L\dAx(t) - -L 1 L* 1 dt 

+ (S - I)dA (t) + dA* (t)L - L* Q S dA (t) - l -L* Q L Q dt 

-iHdt)U(t), (51) 

with initial condition U(0) =1.1 The operators S, L, H , So, L\ and Lq are system operators, where 

• S, L, H, Sq are even (and thus also their adjoints), and 

• L\ and Lq are odd (and thus also their adjoints). 



The operator H is called the Hamiltonian, and it describes the behavior of the system in the absence of 
field coupling. The operators S, L, So, L\ and Lq describe how the field channels couple to the system 
(S and So are required to be unitary). Note that terms involving the creation and annihilation operators 
in (51) ensure a total energy conserving exchange of energy between the system and the field channels; 
for example, an electron may transfer from the field to a quantum dot, and vice versa. Consequences of 
the specified parity of the above operators and the fact that U(0) = I is even is that U(t) is even and 
hence commutes with all the Ito differentials, and, by the quantum Ito rule, is a unitary process (we have 



dA* Lo = -LodA* Q , dA* {t)L x = —LidA*(t), and dB* L = LdB* , see Appendix [AT2| equations ( pPf and 
p4|). 



3.3.1 Heisenberg Picture Dynamics 

A system operator X at time t is given in the Heisenberg picture by X(t) = jt(X) = U(t)* XU(t) and it 
follows from the quantum Ito calculus and the commutation and anticommutation relations arising from 
the chosen parity that 

dj t (X) = Jt (S*XS - X)dA(t) + dB(ty Jt (S*[X,L}) + j t ([L* ,X]S)dB(t) + Jt (£(X))dt 
+dA 1 (t)* Jt (T(X)L 1 - L 1 X)+j t {L* 1 T{X) - XLDdA^+jtiCKX^dt 
+3t(S* XSo - X)dA {t) + dA (t)*j t (S* (T(X)L - L X)) + j t ((L* r(X) - XL*)S )dA (t) 
+j t {£l(X))dt-ij t {[X,H])dt, (52) 

where 

C(X) = L*XL-^XL*L-^L*LX, (53) 
Cl(X) = LxriX^l-^XLxLl-^LxLlX, (54) 
Cl(X) = L*t(X)Lo-^XL*L -^L*L X, (55) 
and in the case of even operators we shall just write £i(X) = L*XLi — \XL\Li — ^L*LiX, (i = 0, 1). 
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The boson and fermion output fields are denned by 

B out {t) = U*(t)B(t)U(t), A out (t) = U*(t)A(t)U(t), (56) 

Ai, out {t) = U*(t)Ai(t)U(t), (57) 

A ,out(t) = U*(t)A (t)U(t), A 0!OUt (t) = U*(t)A (t)U(t) (58) 
and satisfy the corresponding quantum stochastic differential equations (QSDEs) 

dB out {t) = Jt {L)dt + j t (S)dB{t), (59) 

dA out (t) = Jt (L*L)dt + dB*(t)j t (S*L)+j t (L*S)dB(t) + dA(t), (60) 

dA hout (t) = jtiLDdt + dA^t), (61) 

dA 0tOUt {t) = j t (L )dt + jt(So)dA (t), (62) 

dA , out (t) = j t (L*L a )dt + dA*(t)j t (S* L )+jt{L* So)dA a (t)+dA a {t). (63) 

3.3.2 The State 

We define a state E[-] on the von Neumann algebra of observables to be an expectation, that is, a linear 
positive normalized map from the observables to the complex numbers; positive meaning that E[X* JT] > 
for any observable X and normalized meaning E[J] = 1, where / is the identity operator. For technical 
reasons we require the state to be continuous in the normal topology, see for instance [30] . We shall assume 
that the state is a product state with respect to the system-environment decomposition: E[X <g> F] = 
(X)s (F) e, for system observable X and environment observable F. In particular we take (-)e to be the 



mean zero gaussian state with covariance (46) and the choice of N = 1 (the Fermi vacuum). 
We say that the state is even if we have 

E o r = E, (64) 

where r is the parity operator that was introduced in Section |3.3| Specifically, this forces all odd ob- 
servables to have mean zero. In quantum theory, the observable quantities must be self-adjoint operators, 
however, it is not necessarily true that all self-adjoint operators are observable as there may exists so-called 
superselection sectors. In the present case, only the even self-adjoint operators are observables. We need 
to ignore states which lead to unphysical correlations between component systems, this is referred to a a 
superselection principle in the quantum physics literature |26) . We need therefore to restrict our interest 
to even states only. More specifically, we shall assume that the factor states ( • )s and (-)e are separately 
even on the system and environment observables respectively. 

The expected values of system operators X evolve in time as follows. Define 

tit(X)=E[j t (X)}. (65) 



Then by taking expectations of ( 52 ) we find that for even observables X 



fH(X)=n t (C(X)+C 1 (X)+C {X)), (66) 



which is called the master equation, and corresponds to the Kolmogorov equation (13). This may be 
expressed in Schrodinger form using the density operator p(t) defined by fJ>t{X) = tr[p(t)X], which exists 

by our assumption of normal continuity of the state. The density operator is then an even positive 
trace-class operator, normalized so that tr[,o(i)] = 1, satisfying the equation 

p(t)=£,*(p(t))+ J C* 1 (p(t))+ J C* (Pm (67) 

where 

C*(p) = LpL*-^L*Lp-~pL*L, (68) 

Cl(p) - LlpLr-^Llp-^pLrLX, (69) 

£q{p) = L pL* - ^L*L p- ^pL* L . (70) 
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4 Fermion Filter 



In this section we suppose that electrons in fermion channel 0, after interaction with the system, can be 
continuously counted; that is, the observables Ao, ut(s), < s < t, are measured. The problem is, given an 
even state E as outlined above, to determine estimates X(t) of system operators X given the measurement 
record. This is a filtering problem involving a signal derived from a fermion field. As mentioned above 
only the even operators may be observable, and in fact the expectation and conditional expectation of all 
odd operators must vanish identically. 

Mathematically, we wish to determine equations for the quantum conditional expectations 

X(t)=TT t (X)=E[j t (X)\& t ). (71) 

Here, X is a system operator, is the algebra generated by the operators Ao jtm t(s), < s < t, a 
commutative von Neumann algebra, and 7r t is the conditional state. In quantum mechanics, conditional 
expectations are not always well defined due to the general lack of commutativity. However, the conditional 



expectation (71 1 is well defined because jt{X) commutes with all operators in the algebra 'Wf This is called 



the non-demolition property, and is a consequence of the system-field model, where fermion field channel 



serves as a probe, see [7] . The quantum conditional expectation ( 71 1 is characterized by the requirement 
that 

E[j t {X)Z] = E[-ir t (X)Z] for all Z € & t . (72) 



Theorem 1 The quantum filter for the conditional expectation (71) is given by itt(X) — for odd observ- 
ables, while for even observables satisfies the equation 

dir t (X) = TT t {-i[X,H}+C(X) + C 1 {X) + C (X))dt 
^t{L* XL 



n (X)\dW(t) (73) 

where W(t) is a ^ martingale (innovations process) given by 

dW{t) = dY{t) - -K t {LlL )dt, W(0) = 0. (74) 

Proof. We derive the filtering equation using the characteristic function method [35] , U, whereby 
we postulate that the filter has the form 

dn t (X) = T t (X)dt + H t (X)dY(t), (75) 

where Tt and T-Lt are to be determined. 

Let / be square integrable, and define a process c/ by 

dc f (t) = f(t)c f (t)dY(t), c/(0) = l. (76) 



Then cj(t) is adapted to and the requirement (72 1 implies that 



E[X(t)c f (t)]=E[ n (X)cf(M (77) 
holds for all /. We will use this relation to find Ft and Wf 



The differential of the LHS of (77 1 is, using the quantum stochastic differential equation (52) and the 
quantum Ito rule, 

dE[X(t)c f (t)} - E[(djt{X))c f (t)+j t (X)dcf{t) + djt(X)dc f {t)} (78) 
= E[c f (t)j t (C(X) + d(X) + £ Q (X))dt + 3t (X)f(t)c f (t)dY(t) 

+f(t)c f (t) Jt (L*T(X)L - XL*L )dt] (79) 
= E[c f (t)j t (£(X) + d(X) + C (X))dt + f(t)c f (t) 3t (L*T(X)L Q )dt} (80) 
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Now using the property (77) we find that 
d 



dt 



E[X(t)c f (t)} = E[c f (t) n (C(X) + d(X) + Cq(X)) + f{t)c f (t)n(L* T(X)L )) 



Next, the differential of the RHS of (77), using the ansatz (75) and the quantum Ito rule, is 
d 



dt 



E[n t (X)c f (t)] 



E[c f (t)(T t (X) + g t (X)ir t (L*L )) 
+f(t)c f (t)(n t (X)ir t (L*L ) + g t (X)n t (L*L ))] 



Equating coefficients of c/ and fcj in equations (81) and (82) gives the equations 

Tr t (C(X) + d(X) + Co(X)) = T t (X) + g t (X)ir t (L*L ) 

n (L* r(X)L ) = n t (X)7r t (L*L ) + g t (X)n t (L*L ) 

from which the filter coefficients are readily determined, and we deduce the full filter equations 



dn t (X) = n(-i[X, H] + C(X) + d{X) + C Q (X))dt + 
Taking X to be odd and even in turn yields the result. ■ 



7T t (L* T(X)L ) 



(81) 



(82) 



(83) 
(84) 



Tr t (X) \dW(t). (85) 



Corollary 2 Let po be the initial even density matrix for the system, then in the Schrddinger picture we 
may define the conditional density operator p(t) by 7r t (X) = tr[/5(t)JC], and obtain the filtering equation 



d P {t) = (£*(p(t)) + £l(p(t)) + Q(P(t))dt + ( f; K 2frM - P^)) dW ^ 

ti{L Q p{t)L*) 



(86) 



with p(0) = po- 



5 Examples 

In this section we consider several examples drawn from the literature. These examples are special cases 
of the model described above (Figure [5| . 

5.1 Quantum Dot 

We consider a quantum dot arrangement discussed in [3T1 sec. 3], as shown in Figure [2] The left ohmic 
contact L is assumed to be a perfect emitter, which we describe by a fermion field channel A\{i) (for which 
we have dAl(t)dAi(t) — dt), while the right ohmic contact R is assumed to be a perfect absorber, given by 
a fermion field channel A (t) (dA (t)dA^(t) = dt). Current flows through the quantum dot by tunneling. 
The quantum dot is described by annihilation and creation operators c and c* , respectively, satisfying 
{c,c*} = 1 and c 2 = 0. The dot is coupled to the two fermi channels via the operators Li = i^/Zc, 
Lq = iy^fRC, and So = /. Here, jl and 7^ are the tunneling rates across the left and right barriers, 
respectively. We take H = 0, and there is no boson field channel in this example. The parity is defined 
such that c and c* are odd. 



The dynamics of the quantum dot are described by a differential equation for c; from ( 52 ) we have 



dc(t) = - lL + lR c (t)dt + i^EdA^t) + iyfmdAo(t). 



Equation (87) is a linear quantum stochastic differential equation, with solution 



c{t) 



~IL+1R . 



c+ / e" = R( - t - s) (iy^IdA 1 {s)+iy^dA (s)). 
Jo 



(87) 



(88) 
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Note, however, that this system is not Gaussian. The influence of the two fermion fields on the dot can 
be seen in these equations through the stochastic integral terms. 

The output fields are given by 



dAi t0ut (t) 
dA 0>out (t) 



-iy/TL c* (t)dt + dAx(t) 
iy/jRc(t)dt + dA (t) 



and 



dA ,out(t) = i R n(t)dt + i^dA*(t)c(t) - i^c*{t)dA (t) + dA (t), 



(89) 
(90) 



(91) 



where n(t) = c*(t)c(t) is the number operator (an even operator) for the quantum dot. The output field 
components exhibit contributions for the dot and the input fields. Using (52), the number operator n(t) 
satisfies the equation 

dn{t) = ^ L {l-n{t))dt-^ R n{t)dt+iy/^{c*{t)dAi{t)-dA\^ (92) 

We turn now to the expected behavior of the quantum dot system, |37j , |31j . The differential equation 
for the unconditional density operator p{t) is, from (67): 



1 



1 



p{t) = lL { c *p{t)c--cc*p{t) - -p(t)cc*)+7R(cp(t)c* - -c*cp(t) - -p(t)c*c). 



(93) 



The expected number of fermions in the dot is defined by n(t) 
equation 

dn(t) 



dt 



7l(1 - n(t)) - "fRn(t). 



2 rw 2' 
E[n(t)], and satisfies the differential 

(94) 



In steady state, the average number of fermions in the dot is n ss = {lL +7_r), reflecting an equilibrium 
balance of electron flow through the source and sink channels. 

Now suppose that the current in the right contact is continuously monitored; this corresponds to 
the output field observable A-Q tOUt {i). We may condition on this information to obtain an estimate of 
the quantum dot occupation number, n(t) = E[n(t)\^%] — ^t(n). Using (73), the stochastic differential 
equation for n{t) is 



dn{t) = 7i(l - n(t))dt - j R n(t)dt - h(t)(dY(t) - j R n(t)dt). 



(95) 



It is worth comparing the form of this filtering equation to the classical Kalman filter ( 24 ) , and the quantum 
filter for an atom monitored by a boson field, ([T]). While the details of the dynamics differs in these cases, 
the filters share the same structure, with an additive correction term involving an innovations process. 
The 'gain' in this correction term is not deterministic, unlike the case of the Kalman filter (24) for the 
conditional mean £(i). 



5.2 Photodetection 

A photodetector is a sensing device that produces an electronic current flow in response to light incident 
upon it. At the quantum level, a discrete output results from the arrival of a photon. A simple model for 
a photodetector involving both boson and fermion fields is described in [16, sec. 8.5]. In this section we 
use this model to describe the detection of photons scattered from an atom, and then we show how to use 
the information in the electron flow to estimate atomic variables using a fermion filter. 

A schematic representation of the detection of the photons emitted by an atom is shown in Figure [3j 
This figure illustrates that the output boson channel for the atom is fed into the input boson channel of 



the detector. The atom is modeled as a two level system on the Hilbert space C (Section 3.1 ) coupled to a 



boson field (Section 3.2). The atom has lowering and raising operators er_ and cr+ = a*_, respectively, and 



our interest is in the atomic observable n = <j + a- counting the quanta in the atom (0 or 1). The detector 
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is modeled as a three-level system with Hilbert space C 3 coupled to boson and fermion field channels. The 
analogs of the lowering and raising operators are the operators 



vjk = \j){k\, j,fc = l,2,3, (96) 

where |1), |2) and |3) denote a basis for C 3 ; thus \j) = &jk\k), as indicated by the arrows in Figure [3j 

The connection of the atom to the detector via the boson channel is an instance of a cascade or 
series connection, [13], [9], [19]. If the time delay between the components is small relative to the other 
timescales involved, then a single markovian model may be used for the combined atom-detector system. 
In this model, the operators 032 and a\ 3 are odd, while <r_, er + , <7i 2 , <Tn, <j 22 and 033 are even. The 
Hamiltonians for both subsystems is taken to be zero. We take H = (a± 2 a+ — 02 if-) (arising from 

the series connection), S = Sq = I and set the coupling operators to be L = i/k<7- + \/7 Lq — ^/to 032, 
L\ = ^/tT 031. The quantum stochastic equations of motion for the atom are 

da-(t) = ~a-(t) + y/K{2n(t)-I)dB(t) (97) 
dn(t) = -Kn(t)dt-y/K(dB*(t)a--(t) + a- + (t)dB(t)) (98) 

and for the detector 

da n (t) = (ja 2 2(t)+jia S3 {t))dt + v ^f{a + (t)a 12 (t) + a 2 i{t)a-(t))dt 

+^{dB*(t)a l2 (t) + v 21 (t)dB(t)) -^(dAl(t)a 31 (t) - a 13 (t)dA 1 (t)) (99) 

da 22 (t) = -(7 + 70)^22 ^(a+(t)a 12 (t) + a 2 i(t)a-(t))dt 

-^{dB*{t)a 12 {t) - <r 2X (t)dB{t)) - ^(dA*(t)a 32 (t) - a 23 (t)dA (t)) (100) 

da 33 (t) = 7 cr 2 2 (t)dt - 71 a 33 (t)dt 

+V^M(t)a 32 (t)+a23W^o(*)) + V7r(^(*) CT 3i(*) + fi3(*)^ 1 (t)) (101) 



da 12 (t) = --( 7 + 7o)ci2(t)d£+ ^\^l(^2 2 {t) - <Ju(t))<j-(t)dt 

+ 2^y( CT 22(t) - a u (jt))dB{t) - ^foa 13 (t)dA a (t) + ^ a 32 {t)dA\{t) (102) 

da 32 (t) = -^(•j + j + j 1 )a 32 (t)dt-^/K^a 31 (t)(7^(t)dt 

-y/l<ra\{t)dB(t) - v/To (a 22 (t) + a 33 {t))dA (t) - y/jl a 12 (t)dAt(t) (103) 
Ti 

da 13 (t) = — -<r 13 (t)dt + y/ivycT23(t)(T-(t)dt 



+Vl cr 23 {t)dB(t) - ^fo<r 12 (t)dA*(t) - VtT (a n (t) + a 33 (t))dAl(t) (104) 
The number operator for the output of fermion channel evolves according to 

dA ,out(t) = 1o<r 22 {t)dt + ^Tb (dA*(t)a 12 (t) + a 21 {t)dA (t)) + dA (t) (105) 
The mean value n{t) = ¥\n(t)\ evolves according to 

h(t) = -Kfi(t), (106) 

and so n(t) — > as t —¥ 00. Thus initial quanta in the atom are lost to the fields. 

To calculate the conditional mean n(t) — E[n(£)|£%], we will make use of the notation 

= <?ftO a (Jp (107) 

where j, k = 1, 2, 3 and a, f3 = — , + and the fact that atomic operators commute with detector operators. 
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The conditional mean h(t) is given by the system of equations 



dn(t) = -Kh(t)dt+( a2 ? + ~^ -h(t))dW{t) 




(108) 



da 22 (t) = -(j + ^o)o- 22 (t)dt-^(a 12+ (t) + a* 12+ (t))dt-a 22 {t)dW{t) 

da 12 +(t) = + 7 + "fo)vi2+{t)dt - y/ivy a 11+ _(t)dt - a 12+ (t)dW(t) 

da 11+ _(t) = - K& 11+ _(t)dt + 7022+- (t)dt + 71CT33+- (t)dt - a 11+ _(t)dW(t) 

da 22 +-(t) = -(K + 7 + 7o)CT22+-(*)dt-o-22+-(t)<W(t) 



(109) 



(110) 

(111) 
(112) 




da 33+ -(t) = -(« + 7i)^33+- {t)dt + 70^22+- ft) + ( <J l 2+ f ? -a 33+ -(t))dW(t) 



(113) 



(114) 



Note that the filter involves estimates of variables associated with the detector. 

6 Conclusion 

In this paper, using the boson and fermion quantum stochastic calculus we have derived the quantum 
filtering equations for a class of open quantum systems (which may be fcrmionic) that are coupled to both 
bosonic and fermionic fields, for the case where the measurement (driving the filter) is that of counting of 
electrons in a fermionic field. For illustration, the filtering equations were calculated for two examples of 
estimating the number of electrons in a quantum dot coupled to an electron source and sink, and that of 
counting the photons emitted by a two level atom via a photodetector which is modelled as a fermionic 
three level system. For both of these examples we find that the resulting set of filtering equations is 
closed (i.e., the quantum filter is completely determined by a finite number of coupled matrix stochastic 
differential equations). 

A Stochastic Calculus 
A.l Classical 

Let w(t) be a Brownian motion (Wiener process). This means that w(t) is an independent increment 
process and w(t) — w(s) is a Gaussian random variable with zero mean and variance t — s. Suppose that 
a(t) is an adapted stochastic process, that is, a(t) is independent of w(r) for all r > t; in particular, the 
Ito increment dw(t) — w(t + dt) — w(t) is independent of a(t). The Ito stochastic integral of a with respect 
to w is defined as a limit involving forward increments: 



where, sq = < si < s 2 < ■ ■ ■ < t. Since the Wiener process has zero mean, so does the stochastic integral: 




(115) 





(116) 



where a and j3 are adapted. Then the Ito product rule says that 




(117) 
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which is the sum of a stochastic integral and conventional (say Lebesgue) integral, called the Ito correction 
term (the last term on the right hand side). The Ito correction term is not present in the usual product 
rule. 

Stochastic integrals and the product rule are often expressed in differential form. Indeed, we may write 

dl = adw and dJ = fidw, (118) 

so that 

d(IJ) = (dl)J + I(dJ) + (dl)(dj) (119) 
= Jadw + If3dw + a/3dt, (120) 

where we see that the correction term arises from the Ito rule 

dw(t)dw{t) = dt (121) 

which is very useful in calculations. 



A. 2 Quantum 



Let B(t), B*(t) be boson annihilation and creation operators, as discussed in Section 3.2 For simplicity 
here we ignore the number (counting) field operators and assume M = and N = (vacuum field state). 

Let a>i(t) and a 2 (t) be operator-valued adapted processes, i.e. independent of future field operators. 
Quantum Ito stochastic integrals are defined in terms of forward increments: 

I(t) = [ ax(s)dB{s) + [ a 2 {s)dB*{s) 
Jo Jo 

« Y,*i(sj)(B(s J+1 )-B(s 3 )) + Y,a2(s 3 )(B*(s ]+1 )-B*( Sj )). (122) 

j 3 

The expected value of the above stochastic integral is zero in the vacuum state \(j>): K,p[I(t)} = 0. 
The quantum Ito rule is expressed in terms of four products 

dB(t)dB(t) =0, dB(t)dB*(t) = dt, dB*(t)dB(t) = 0, dB*(t)dB*(t) = 0. (123) 

This Ito table is valid for vacuum and coherent field states. Ito tables for squeezed and thermal field states 
have more non-zero terms, as mentioned in Section 3.2 An important property is that for an adapted 
process a(t), we have 

[a(t),dB(t)] = = [a(t),dB*{t)]. (124) 
Now suppose we have two quantum stochastic integrals, I(t) defined above and 

J(t)= f /3 1 ( S )dB(s)+ f p 2 {s)dB*{s). (125) 
Jo Jo 

The product rule is 

d(IJ) = {dI)J + I{dJ) + (dI){d.J) (126) 
= aiJdB + a 2 JdB* + I(3idB + Ip 2 dB* + aifcdt (127) 
= (aiJ + I0x)dB + (at l J + iP2)dB* + a 1 fodt; (128) 

that is, 

I(t)J{t) = f (ai(s)J(s) + I(s)^(s))dB{s) + J (a 2 (s)J(s) + I(s)/3 2 (s))dB*(s) 

a 1 (s)f3 2 (s)ds. (129) 
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The last term in this expression is the Ito correction terms, and arises from the non-zero Ito product 
dB{t)dB* (t) = dt. Note that the order is important in these expressions, since the expressions involve 
quantities that need not commute. 

Quantum stochastic integrals with respect to fermion fields may also be defined, [2], [Hj- However, 
some extra effort is required to keep track of antisymmetric tensor products and parity, matters that are 
explained in Appendix |B| Here we mention that for an odd adapted process j3{t) and a fermion field A(t), 
A* (t) we have the anti-commutation relations 



B Mixed Boson-Fermion Systems and Generalized Antisymmet- 
ric Tensor Product of Operators 

The laws of quantum physics dictate that fermionic operators from independent systems must anticommute 
with one another. That is, if Xj is any fermionic operator on a system j with Hilbert space Sjj for j = 1,2 
then we must have {X±, X 2 } — X\X 2 + X 2 Xi — 0. Realizing this property is not trivial because typically 
when one has a composite system with Hilbert space Sji ®f) z (here ®, depending on the context, denotes 
the tensor product between two Hilbert spaces or the tensor product between two Hilbert space operators) 
the operator X\, originally defined only on Sji, would be identified with its ampliation X\ ® I on the 
composite Hilbert space, whereas X 2 would be identified with / ® X 2 . This is problematic for fermionic 
systems since then we would have that {Xi®I, I<EiX 2 } = 2X 1 (g)X 2 , which is not necessarily for arbitrary 
fermionic operators Xi and X 2 . To resolve this issue and get these operators to anticommute correctly, 
the Hilbert spaces of fermionic systems are endowed with some additional structure and the usual tensor 
product between operators must be replaced with the antisymmetric tensor product between operators. 
We briefly explain this below, for a more detailed account see, for instance, [2]. 

Let F) be the Hilbert space of a fermionic system. Then it is required that there exists an orthogonal 
decomposition Sj — $) + ©fj_ where the subspace Sj + is called the even subspace, while is called the odd 
subspace. An operator X on f) is said to be even if it leaves the even and odd subspaces invariant XSj± C Sj± 
and odd if it maps vectors of the even subspace into the odd subspace and vice- versa X$)± C Sj T . In 
particular the fermionic operators (or degrees of freedoms) are odd operators. 

Let P+ and P_ denote the orthogonal projection onto the odd and even subspace, respectively. Define 
the linear parity operator 9 = P + — P_ and the linear parity superoperator t(X) = 6X9 for any bounded 
operator X on . It is easily checked from the definition that 9 and r have the following properties: 

1. 9* = 9. 

2. 9*9 = 9 2 =1. 

3. r(XY) = r(X)r(Y) 

4. t(X*) = t(X)* 

5. If X is even then t(X) = X, while if X is odd then t(X) = -X. 

Note that properties 3 to 5 imply that r is an automorphism on £>(f)) (the space of all bounded operators 
on Sj). For an operator X with definite parity (i.e., it is either even or odd) we define the binary- value 
parity functional 6(X) as S(X) = if X is even and S(X) = 1 if X is odd. Thus for any operator with a 
definite parity we have that t(X) = (— l) st - x ^X . In general an operator X has a decomposition into odd 
part and even parts as X = X + + X-, where X + and X_ are even and odd operators, respectively, given 



{0(t),dA(t)} = O = {l3(t),dA*(t)}. 



(130) 



by: 



X- 



P+XP+ + P_IP_ 
P+XP- + P-XP+ 
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Thus 

t(X) = T (X+) + t(X_) = (-1) S ^ X ^X + + (-l) s ( x -1X_ =X+-X_. 

We can now define the composition of independent fermionic systems (f)i,#i) and (^2,^2)- The com- 
posite system has Hilbert space SJ12 = Sji (8)$)2 and parity operator 9 12 — &i ® 82 so that its even and odd 
subspaces are -H12+ = H 1+ ® H 2 + © Hi- (g) H 2 - and H12- = Hi- ® H 2 + © H\ + ® H 2 ~- The key question 
is how to ampliate the action of operators Xj acting on the individual Hilbert spaces Sjj to the composite 
Hilbert space, so as to satisfy the desired fermionic anti-commutation relations. For this we define the 
antisymmetric tensor product © between Xi and X2 as 

Xi®X 2 =X 1 <2> X 2+ + Xi9i © X 2 -. 

The ampliation of an operator X\ of the first system is Xi®I = X\ ®I, and the ampliation of an operator 
X 2 of the second system is I®X 2 = I®X 2 + + 9i®X2-, such that X X ®X 2 = (Xi®I)(I®X 2 ) ■ In particular 
if both Xi and X 2 are odd then {Xi<E)I, I(g)X 2 } = (the ampliations anticommute), while if at least one 
of them is even then [Xi<g>I, I&X2] = (the ampliations commute). Note that the ampliation preserves 
the original parities of the system operators, e.g. Ti2(-f<8>^2) = 9i2(I<8>X 2 )9i2 = I<E)T2(X 2 ). 

Now, if (# 3 , 6*3) is a third fermionic system we can show that (Xi®X2)®X 3 = Xi®(X2<§)X 3 ) so that the 
antisymmetric tensor product is associative and Xi®X 2 ®X 3 is unambigously defined. By repeating the 
above procedure we can define the composition of an arbitrary number of systems (fjj, 9j), j — 1, 2, . . . , n 
such that ®j—iXj is well defined for arbitrary operators Xj on fjj. For simplicity we identify each Xj 
with its ampliation to ©™ =1 "Hj with respect to the antisymmetric tensor product <g>, so that if j 7^ k then 
{Xj,Xk} = whenever both Xj and X^ are odd, while [Xj, Xk] = whenever one of them is even. 

In the setting of this paper, we also consider composite systems consisting of fermionic and non- 
fcrmionic (bosonic) sub-systems. The composition can be described in the framework of the anti-symmetric 
tensor product by endowing the non-fermionic systems with a trivial parity structure. Indeed according 
to quantum physics, bosonic operators from one system commute with all operators (even or not) from 
the other systems, so they can be interpreted as "even" operators. This can be achieved by defining the 
parity operator of a bosonic system as 9 = I (i..e., take P + = I and P_ = 0) so that t(X) = X for 
all bosonic operators. In this way we can extend the definition of ® mutatis mutandis to mixtures of 
fermionic and non-fermionic systems in the way we had done above for fermionic systems. For instance, 
let S)i and Sj 3 be the Hilbert spaces of two non-fermionic systems and $j 2 is the Hilbert space of a fermionic 
system. Consider the composite system on Sj\ ©^2 and let Xj be an operator on Sjj. Then we have that 
Xi®X 2 = Xi ® X 2 + + Xi © X2- = X\ © X 2 - Similary, for the composite system on S32 © -£>3 we have 
X2<8>X3 = X 2 © X$. That is, the generalized antisymmetric tensor product between two operators from 
distinct systems of which one is fermionic and the other is not, reduces to the usual tensor product between 
operators. On the other hand, if the two operators are from two distinct fermionic systems then it reduces 
to the antisymmetric tensor product for operators of fermionic systems. This is exactly how it should 
be: non-fermionic operators from one system commute with all operators from the other systems, and the 
fermionic ones anti-commute with fermionic operators from the other systems. By inspection it is easy 
to see that the generalized antisymmetric tensor product defined in this way has the associative property, 
hence it is unambiguously defined for systems that are the composite of any finite number of fermionic 
and non-fermionic systems. 
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